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Abstract 

To understand the direct J/il^ + Xnon-cc production mechanism in e^e~ annihilation, in this 
work, we propose to measure the inclusive J/^ plus light hadrons (LH) production at B-factories 
and present a detailed study on its QED production due to ip{2S) feed-down, where the ip{2S) 
are produced in e"'"e~ — )• 'ip{2S) + 7 and e^e~ — )• 'ip{2S) + ff, f = lepton,lightquark, and QED 
contribution to direct J/^ -|- qq production with q = u,d,s quark. We find that the QED con- 
tribution is huge in the whole phase space region, but can be reduced largely and is in the same 
order as the QCD contribution when a suitable cut on the angel Oj/^ between J/ip and the e~^e~ 
beam is made. In this way, the cross section of J/V' + LH QCD production can be obtained by 
subtracting the QED contribution from the experimental measurement on inclusive J/ip plus light 
hadrons. To help to remove the QED background, we also calculate the angular and momentum 
distribution of J/^ in the QED contribution. 

PACS numbers: 12.38.Bx, 13.66.Bc, 14.40.Pq 
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I. INTRODUCTION 



The development of the non-relativistic QCD (NRQCD) effective field theory [l| provides 
a powerful tool to study the production and decay of heavy quarkonium states that are con- 
stituted by one heavy quark (Q) and one heavy anti-quark Q. The virtual difference between 
NRQCD and the conventional color-singlet model (CSM) is that it allows the contribution 
of QQ state in the color-octet (CO) configuration at short-distance which finally evolves into 
heavy mesons through emission of soft gluon(s) non-perturbatively. This is refereed as the 
CO mechanism (COM). The role of the COM has been extensively studied in various high 
energy environments, for reviews see Ref. 2|. 

Among them, the J/ if) production in e^e~ annihilation at B-factories (the Babar and 
Belle) have attracted considerable solicitude in recent years. Experimentally, the cross sec- 
tion for inclusive J /ip production was reported by the Babar [3i] and Belle collaborations 
in 2001. Belle collaboration further divided the inclusive J/ip production rate into two 
pieces: (a) e+e" — )■ J jip + cc part^, (b) c'^e' J /ip + Xnon-cc part, and measured them 



separately 



6j. The latest results reported by the Belle are [g 

a(e+e" ^ J/^ + X) = 1.17 ± 0.02 ± O.OTpb, 



(la) 



a(e+e" J/i> + cc + X) = 0.74 ± 0.08^JJ:[!^pb, (lb) 
a(e+e- ^ J /i) + Xnon-cc) = 0.43 ± 0.09 ± 0.09pb. (Ic) 
In the case of Jlib + cc production, where the the CO contribution is found to be very 

n 

small (91], there were large discrepancies between Belle results and NRQCD predictions at 



leading order (LO) in as and v js-lld], where v is the relative velocity between c and c 



in the meson rest frame. These puzzles are now large 
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y resolved after taking into account 



12] and the relativistic corrections 



;he next-to- leading order (NLO) QCD corrections 

1^. In contrast, in the J / + X^on-cc production case, the contribution of the CO e"'"e~ — )■ 
cc{}Sq, '^Pj) + g 14j-ll7l| process is expected to be significant and even larger than that 
of the CS e^e~ — )■ cci^S]) + qq process. At LO, the cross sections of the CO and CS 
processes are predicted to be about 0.3 ~ 0.8pb and 0.2 ~ 0.3pb [8|, 118|], respectively. 



Here the J + cc includes both the exclusive double charmonium production process and the J /ij) 
production process in association with charmed hadrons , and the_cross section for double charmonium 
production process was also reported by Babar collaboration later 
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Recently, the k- factor of their NLO QCD corrections were found to be about 1.3 19|, |20| 
and 1.9 [21] correspondingly, and what's more the relativistic corrections can also enhance 



the LO CS result by about 20% - 30% [221, l23|. Then up to the NLO of a, and t;^, the 



cross section of the CS contribution itself can reach about 440 ~ 560fb [22]. which almost 
saturate the Belle measurement and leave very little room for CO contribution. So in 
e'^e~ J/ip + Xnon-cc process, there is a large conflict between NRQCD prediction and 
Belle current measurements. By setting the CS contribution into zero, the upper limit of 



the CO matrix elements is obtained 



2l| 



(0|C(^5o')|0)^/'^ + 4.0(0|C)(=^P^)|0)-^/^/m^ < (2.0 ± 0.6) x lO-^GeV^^ 



(2) 



However, in some other processes of J/ip production, the recent theoretical calculations 
shown that the CO contribution is important. For example, (a) for J ftp production from 
Z decay, the CS result at QCD NLO 2^, can only account for one-half of the experi- 
mental data and the other half might be attributed to the CO contribution; (b) for J/ip 



production in T decay, there is large gap between CS contribution [25| and the experimental 
result; (c) for J/ip photoproduction at HERA, the transverse momentum (pt) distribution 
and the polarization parameters of J/ip can not be well described by the CS channel at 



QCD NLO alone as well [261, ui, and the NRQCD prediction that includes both the CO 



and CS contributions can give a well description of the J/ip pt distribution when the NLO 



QCD corrections are taken into account 



281 ]: (d) for J/ip hadropro duct ion, despite of the 



huge NLO QCD corrections 29143 1|. the CS contribution still can not explain the experi 



32-|35|. By fitting the J/^p 



mental measurements, and the role of the COM is significant 
hadroproduction data with the complete NRQCD results at QCD NLO, including both the 
CS and CO contributions, two different sets of constraint for the CO matrix elements are 



obtained, which are {0\O{^Sl)\0)'^/^ + 3.9{0\O{^P^)\0)^/^ /m^= 7.4 x lO^^GeV^ y, and 



35j . Although these two results 



{0\O{^S^)\0y/'^ + 3.9{0\O{^Pl)\0y/^ /ml = 2.4x lO-^GeV' 
are not consistent with each other, both of them exceed the upper limit given in Eq.([2]). In 
particular, the former one is three times larger than the limit in Eq.(l2]). These studies yield 
almost completely opposite conclusion about how large the CO contribution is, or in other 
words, how large the values of the CO matrix elements could be. 

After comparing the results of the Babar and Belle with the theoretical calculation metic- 
ulously, we find that there are some uncertainties which can potentially have large impact 
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on the current conclusion. One is that the Babar measurement on J lib inclusive production 

n n 

cross section is about two times larger than that of Belle [6]. If we subtract the mea- 
surement of the Belle (j{J /ip + cc + X) = 0.74pb, which is well understood theoretically, 
from Babar result, there will be enough room left for CO contribution. One possible reason 
for the different results of the Babar and Belle is that they use different methods to select 

n 

the data. Another uncertainty is that, in the latest measurement of the Belle [6[, they only 
select the event that includes at least five charge tracks in the final states, and make no 
corrections. This means that all events that include zero or two charged light hadrons, such 
as J /ip + m{'K^'K^) +ra7r° for (m = 0, 1; n = 0, 1, 2 . . .), are excluded. From the point view of 
quark-hadron duality. Belles measurements do not include the whole NRQCD predictions. 
It may has little influence on the measurement of a{e^e^ — ?■ J /ip + cc) [6], but large influence 
on that of cr(e+e~ J /ip + Xnon-cc) from the non-perturbative hadronization mechanism 
of gluons. To reduce the uncertainties mentioned above and understand the J /^J + Xnon-cc 
production mechanism, we suggest to measure the cross section of J/'?/' + light hadrons (LH) 
production by the Belle and Babar collaborations with the same kinematic criteria, which 
can be compared with the theoretical prediction directly. 

Besides the interesting conventional QCD contribution, there are also large QED back- 
grounds due to ipi^S) —7- J jip + vrvr^, where V^(2S') is produced in the initial state radiation 
(ISR) process e^e~ — j- ip{2S) + 7 and higher order QED processes e+e^ — )■ ip{2S) + // (/ 
can be lepton or light quark), and direct J/ip production in the e'^e' — t- 27* — )■ J/ip + qq 
process with q = u,d,s quark. To help to remove them, in this work, we will present a 
detailed study about the ipi^S) and J/ip + qq productions in the QED processes and their 
influence on the J /ip + LH measurement. 



II. FRAMEWORK OF CALCULATION 



For the process of J /ip + t^ti production from ip{2S) feed-down, the Feynman amplitude 
Ai can be generally written as: 

^ V'(25') can also decay into J/ip + rj. However the branching ratio is more 15 times smaher than the 27r 
channel, so we do not take it into account in our calculation. 
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where A4t.^'^^\P2s) and _A/fW'(^'^)~^'^/'^+'^'^) are the Feynman amphtudes for il){2S) production 
with momentum P2S and ip{2S) decay into J/ip + rnr respectively, and F is the total decay 
width of ip{2S). Using narrow width approximation 

li,„ I = -'^P's - Ml,) 

it is straightforward to obtain the expression for the corresponding cross section which is 
factorized as the product of the cross section of ipi^S) production and the branching function 
of ^{2S) J/^ + tttt: 



a 

8s 



where LIPSi is the phase space of ipi^S) production, LIPS2 is the phase space of ip{2S) 
decay into J /ip + im, and J = 1 is the spin of ip{2S). 



We use the effective Lagrangian that is constructed in Ref. [39l| to describe ip{2S) — )■ 
J/ip + nn. The amplitude J\A^'^^'^^'^^"^^^~^'^'^^ can be read directly from the Lagrangian 

where m^^ = [pi + ^2)^? is the mass of vr meson, and v = (1, 0) in the rest frame of 
ip{2S). In their convention, the vr decay constant Fq ^ 93MeV. The coupling constant 
(72 ~ 0, because it is strongly suppressed by the chiral symmetry breaking scale over m^. By 
fitting the distributions of 771^^^^ and cos 6'*, which is the angel between J/ip and Tr"*" in the 
rest frame of iIj{2S), in the decay of il){2S) J / ip + t^^ , the BES Collaboration obtained 



two set results for — and — 
9 a 

they then obtained that^ 



40j. Together with Br(V'(25) ^ J/^p + tt+tt-) = 33.6% |4l|. 



0.322, — = -0.49, — = 0.54, (7) 
9 9 



or 



g = 0.319, ^ = -0.347, = 0. (8) 

For the processes considered, is common, so we essentially only need to 

compute Ai"^^ . In the non-relativistic limit, for the QED process of e+e^ — ?• ^(25*) + X the 



These parameters can also well reproduce the decay width of ip{2S) — >■ J/ip + tt^tt". 
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AAA/W: 



-* KA/WWWWV7 

(a) 

ij{2S) 

-> fVAAAAjcr^ 





Jli> 





FIG. 1: The typical diagrams describing the 'ip{2S) feed-down background. 



factorization formula in CSM and NRQCD are equivalent, and the amplitude M^^"^^^ can 
be written as: 

= ^ J2 ^^)(3^; 3j\l)M{e^e- ^ c,i^, s,) + c,{^, s^) +X) (9) 

where Ai is the standard Feynman amplitude for e~^e~ — )■ Cj(^,Si) + 0^(^,52) + X, 
(3z;3j|l) = I/a/ZVc and (si;s2|l 5*2) are the SU(3)-color and SU(2)-spin Clebsch-Gordan 
coefficients for cc projecting on the CS spin-triplet S-wave state. The projection of Dirac 
spinors can be re-expressed as: 

1 



'iS,)iP2S + M2s). 



(10) 



C25 can be related to the ipi^S) wave function at origin by C25 = ^|-R2S'(0)|. And |i?25(0)| 
can be obtained from potential model calculation or can be determined from ip{2S) decay 
into e^e~ with 

16a2|i?25(^"' 



T{ip{2S) e+e" 



(11) 



III. THE FEED-DOWN BACKGROUND FROM e+e ^{2S) + 7 

The typical Feynman diagrams for the ISR process e"'"e^ — )■ ip{2S) + 7 followed by 
ip{2S) — 7- J/ip + TTTT are shown in Fig. (la). Using the formula introduced in Eq.(5-8), we 
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compute \Ai{e~^e — )■ ip{2S) + 7)^ analytically and obtain 

sr (-1 + r + (-1 + r) (-1 + 4re) x^) ^ 
-8re + 16rre + 32r2-4 (r + 2re) (-l+4re) + (-1 + r)^ (1 - 4re)^ a;^) (12) 

where Cc = |, r = re = X2 = cos(6'^(25)) and 6*^(25) is the angel between il){2S) 
and the e^e~ beam. In the limit of r^. = 0, Eq. f|T2l) can be simplified as: 

Setting M2S = 3.686GeV, me = O.SlMeV, a = and using 1(^(25) -> e+e") = 
4.30keV, we get 

(j(e+e- ^ ^(2^) + 7) = 13.22pb. (14) 
And the feed-down production 

cr(e+e^ ^ ^'(25) + 7) x 5r(^/'(25) ^ J/^ + tttt) = 6.79pb, (15) 

which, as expected, is huge. This is because in the limit of mg — )■ 0, there will be coUinear 
singularities in |A^(e+e^ — )■ ip{2S) + 7)^ in Eq. (fT3|) at X2 = ±1 points. The angular 
distribution ^'^^^'^'^ 'dxf^^^^'^^ shown in Fig. [2]. It can be found from Fig. [2] that the 
differential cross section drops down very fast when ip{2S) goes off the beam line a little. 
If we make a cut on X2, i.e the angle 6^^(25), the cross section will be reduced largely. The 
cross sections in different cut conditions are given below: 

a(e+e- ^ ^"(2^) + 7) x Br{ij{2S) J/^ + tttt) ^ = 1.48pb, (16a) 

(T(e+e~ ^ ^(25) + 7) x Br{tp{2S) -> J/'ip + tttt) ^ = 0.99pb, (16b) 



cj(e+e- ^ ip{2S) + 7) x Br{ij{2S) J/ip + tttt) ^ = 0.71pb. (16c) 

6<^V(2S)<-6- 

Let p^^^ denotes the four-momentum of J /ip in the rest frame of ip{2S), then |pjy^|/i?}y^, 
the three- velocity of J/V', ranges from to 0.15, which is much smaller than that of '4'{2S) 
in the center of mass frame (CMF) of e^e~ collision, which is about 0.78. So the angular 
distribution of J /ip can be obtained approximately by setting 6^(2S) = ^j/ip, where 9j/^ is the 
angel between J/ip and the e~^e~ beam. However, such an approximation may not be good 
enough here, because the cross section of J/ip + tttt produced from the feed-down of 4'{2S) 



FIG. 2: The angular distribution of il){2S) in the ISR process e+e i^{2S) + 7, where X2 
cos(0^(25))) ^iid ^^,(25) is the angular between Tp{2S) and the e+e^ beam. 



ISR process (Eq.lfTSl)) is more than 10 times larger than that of the CS QCD process |22|], and 
a tiny difference may potentially result in a considerable effect. In this work, we calculated 
it directly. In the CMF, X2 = cos{9j/^), can be expressed as: 



(17) 



|ki||Pj/^l 

where ki is the four- momentum of e"*", L is the Lorentz transformation from 'ip{2s) rest frame 
to the CMF. p'j^^ is the J/ip four-momentum in the CMF. To do the calculation, the formula 
for the decay ip{2s) — ?■ J/ip + vrvr and ip{2s) production are placed in the numerical phase 
space integration program generated by using the Feynman Diagram Calculation (FDC) 
package 42], in which the Lorentz transformation and the cut conditions are employed in 
the numerical calculation. The decay ip{2s) J / ip + rni is calculated by using Eq.(|6]) with 
= M^- = 140MeV and M^o = 135MeV. When parameter set in Eq.([7D are used, the 
cross sections in different cut conditions are 



cr(e+e ^ (J/^ + 7r7r)^(25) +7) 

a(e+e" {J /i) + 7r7r)^(25) + 7) 
cr(e+e" -^{J/il) + 7171)^23) + 7) 



I8<6'j/^<1^ 



f <^'j/v.<-r 



s<»//v.<— 



= l.Slpb; 

0.99pb; 
0.71pb. 



(18a) 

(18b) 
(18c) 
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Alternatively, if we choose the parameter set in Eq.(|8]), the corresponding cross sections 
become: 



a(e+e- ^ ( J/V^ + vr7r)^(25) + 7) , = l-52pb; (19a) 

fs<S.j/i,<^ 



(T(e+e- ^ ( W + 7r7r)^(25) + 7) « = l.OOpb; (19b) 



f <^'J/V)<— 



= 0.71pb. (19c) 



The numerical results in Eq. fll6|18|[T9|) show that for J/ip + tttt production from the ISR 
'ip{2S) feed-down process the approximation 

da{e+e- {J/^j + 7r7r)^(25) + 7) _ da{e+e- ^{2S) + 7) x Br{ip{2S) J/^ + tttt) 
dcos{6j/^) d cos{6^(^2S)) 

(20) 

holds very well in the range of 7r/9 < 9j/^ < 8tt/9 at y/s = 10.6 GeV, and the J/ip 
angular distribution is almost not dependent on the details about how 4'{2S) decays into 
J/ip + Tiii. Hence the angular distribution of J/ip can be safely obtained by using the angular 
distribution of ip{2S) in the interval 7r/9 < Ojj^ < 87r/9 with an additional renormalization 
factor of branching ratio of ip{'2,S) J /iIj + tttt. 

Because the energy difference between ip{2S) and J /ip is at the same order as the energy 
of the soft gluon emitted from the CO cci^Pj,^ Sq) states [l|, which is of mcv"^ order, there is 
a large overlap between the kinematic region of the J/tp coming from ISR ip{2S) feed-down 
and that of the J/ip produced in the CO process. To measure the CO J/ip production, it 
is helpful to know the momentum distribution of J/ip production in the feed-down from 
the ISR ip{2S) process. We calculate it numerically with different cut conditions of Oj/,^ 
by using the two set of parameters in Eq. flTfHj) . and the results are given in Fig. [3(a) 3(d)] . 
The results in Fig. [3(a) 3(d)| show that similar to the angular distribution, the momentum 
spectra give almost same results for those two different parameter sets. Hence for simplicity, 
we only choose the parameter set in Eq.(I7]) in the following calculations. 



IV. BACKGROUND FROM HIGHER QED PROCESSES 



Dackground coming from higher QED processes e^e — )• 



In direct J/ip production, the 
J/i^ + f + f is also considerable 
the feed-down background from e~^e~ — )■ i'i^S) + f + f can not be ignored too. The typical 



38| . where / can be lepton or light quark, and therefore 
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2,5-1 

12 - 




Pj/„ (GeV) pj,^ (GeV) 

(c) (d) 

FIG. 3: The momentum spectra of J /tp produced from the feed-down of ISR ip{2S) process in 
different cut condition of 6j/^ by using two different sets of parameters in Eq.([7]) (sohd hne) and 
Eq-dH]) (dashed hne) to describe ■ip{2S) J ftp + vrvr. 

Feynman diagrams for f ^ e are shown in Fig.(lb) and Fig. (Id). When / = e, there are 
additional t-channel diagrams, the typical one of which is shown in Fig.(lc). Because of this 
t-channel enhancement, the cross section for / = e is expected to be much larger than f ^ e 
case. We will discuss / = e and / 7^ e cases separately in the subsections. At this order, in 
addition to the ip{2S) feed-down, there is also sizable QED contribution from direct J/ip + qq 
production with q = u,d,s quark, about which we will discuss in subsection C. 
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A. The Feed-Down Background From e+e V'(25') + e+e 

According to the interaction type of the initial e+e~, we divide the e+e~ — >■ ■0(2-5') + 
e+e~" process into three part: the t-channel part (Fig.(lc)), the two-photon channel part 
(Fig. (lb)), and the s-channel part (Fig. (Id)). It is easy to check that the Feynman amplitude 
for each part itself is gauge invariant. Compared to cross section a'^ for the t-channel part, 



the cross sections for the two-photon part a"^^ and the s-channel part are suppressed by 
the factors — and '^'^^'^^ ln~^(^^) respectively, which are about 10~^ and 10~^ orders 
accordingly at ^/s — 10.6GeV. Choosing the same values for the parameters as in the ISR 
process, we obtained 

a^{e^e~ [Jj-^ + 7r7r)v,(25) + e+e") = 0.50 pb; (21a) 

(7^P(e+e- ^ (J/V' + 7r7r)v,(25) + e+e") = 4.8 x lO"' pb; (21b) 

(T^(e+e- ^ ( J/V' + 7r7r)^(25) + e+e") = 8.5 x 10"^ pb. (21c) 

which are consistent with the qualitative estimation. The contribution of the s-channel part 
is only ~ 1 fb order, which is about three times order less that the t-channel contribution, 
so we drop it in the later analysis. 

If we make the same cut on the ^j/^, cr^ and a^^ both drop down largely too: 



^T(TP)^g-hg- ^ ^jj^ ^ 7r7r)^(25) + e+e") ^ = 0.11(0.019)pb; (22a) 

^T(TP)^g+g- ^ ^ji^ ^ 7r7r)v,(25) + e+e") = 0.059(0.013)pb; (22b) 



\<Q.JIih< 



J/V^-g- 



crT(TP)(g+g- ^ ( j/^ ^ 7r7r)^(25) + e+e") ^ = 0.039(0.010)pb. (22c) 

The angular distribution of ip{2S) in the t- and two-photon channel parts are shown in 
Fig.[4]. We find that in 7r/9 < 6j/^ < Svr/Q region the angular distribution of the J/ip 
production from feed-down can be obtained by using that of ip{2S) as well with an additional 
renormahzation factor of Br{ip{2S) — >■ J/t/j + nn). Unlike the ISR process, pj/^ ranges from 
to 4.7 GeV, and the momentum spectra of J/t/j for the the t- and two- photon channel 
parts are shown in Fig. [5]. 

We also calculate the interference between the t-channel part and the two-photon part 
and find it is very small. The cross section of the interference part for {J/iP+t:t:)^(^2S)+^~^^~) 
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(a) (b) 

FIG. 4: The angular distribution of tp{2S) produced through t-channel (a), and two-photon 
channel (b) in the e+e" — >• ip{2S) + e^e~ process, where X2 = cos(0^(2S))) and 0^(25) is the angular 
between 'ijj{2S) and the e"*"e~ beam. 



is about — 20fb in the whole phase space region. After including the interference part, the 
total cross section in different cut conditions of 6{J/tp) are 



a{e^e ^ ( J/z/^ + 7r7r)^(25) + e+e ) 
(T(e+e" {J/i{j + 7171)^(^23) + e+e") 



ig<yj/v.<T8- 



g <6'j/^<-g- 



= 0.12pb; 

0.070pb; 
0.047pb. 



(23a) 

(23b) 
(23c) 



a(e+e ^ (J/^A + vr7r)^(25) + e+e ) 

The angular distribution of ip{2S) and momentum distribution of J/ip for the whole pro- 
cess can be approximately obtained by adding the t- and two-photon channel contribution 
together respectively, because the interference effect is very small. 



B. The Feed-Down Background From e+e —?- ij{2S) + ff {f ^ e) 



The process e~^e — )■ ipi^S) + f f (f e) has been fully studied in Ref . |43| . We also 
compute it independently and obtain consistent results: 



J2 a{e+e- ^ ( J/^ + 7r7r)^(2s) + //) = 0.026pb. 



(24) 
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0.08- 

> 

S 0.07 ■ 
g 

■S 0.06 ■ 



0.01 - 



(a) (b) 

FIG. 5: The momentum distribution of J/^[J produced through the t-channel (a) and two-photon 
channel(b) in the e"'"e~ — >• ip{2S) + e~^e~ process. 



If we make the same cut on the 6j/^, the cross section becomes: 



f=fl,T,U,d,S 

J2 f^(e+e- ^ (J/^ + 7r7r)^(2s) + //) 

f=fl,T,U,d,S 



Tg<fj/V< — 



^<Sj/4,<^ 



(J/V' + vr7r)^(2s) + //) 



= 0.020pb; 

O.OlSpb; 
O.Ollpb. 



(25a) 

(25b) 
(25c) 



f=fl,T,U,d,S 

The cross sections in different cut regions are within only about 0.020pb, which are about 
4 ~ 6 times less than those in the e~^e~ — )■ ip{2S) + e'^e~ process, so small that we will not 
present further analysis here, and recommend Ref. 43|] for more detailed results. 



The Background From e^e 



The Feynman diagrams for the e^e^ — ?■ J/tp + qq process are similar to those for e^e~ — )■ 
ipi'^S) + qq process. Since in the V'(2S') production process the contribution of the s-channel 
diagrams can be ignored, for the same reason, we will not consider it here too. The cross 
section of the e+e" — )■ J/ip + qq has been calculated in Ref. |43l| . which is also considerable. 
Using the method introduced in Ref|43|, we calculate the cross section with different cut 
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(a) (b) 

FIG. 6: The angular (a) and momentum (b) distributions of J /tp in the process of e'^e" 
J/V' + qq-, where x'2 = cos(0j/^), and Ojj^ is the angular between Ojj^ and the e^e~ beam. 



conditions of 9 



q=u,d,s 



q=u,d,s 



7r ^/) ^ 177r 
lg<y.J/^<l^ 



■9<yj/^< — 



-- 0.071pb; 

0.052pb; 
0.039pb. 
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q=u,d,s 

The J/ip angular and momentum distribution are shown in Fig. [6]. Note the difference 



between our results and those in Ref. 



43| is due to the different choice of the parameters and 



the amount of data samples used in the i?- value curve 41 1 



V. DISCUSSIONS AND CONCLUSIONS 



Summing up the feed-down contribution from the ISR and // processes and the contri- 
bution of direct J/ip + qq production, the total QED background cross section are about 



(jQEDie^e- J/i) + LH) = 7.46pb, 



(27) 



which is more than one order of magnitude 
QCD production e+e" ^ J/^j + LH [sI, Q- 



arger than the cross section of the conventional 
2l|. Such huge background make it difficult to 
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measure the QCD contribution in the whole phase space region. However, the background 
in the off beam region will drop down deeply. The cross section of the QED background in 
different cut regions are: 



(^QEoie^e ^ J/i) + LH) 

(^QED{e^e' J/i) + LH) 
(TQEoie^e- J/ip + LH) 



T8<»J/V< — 



9 <6'j/V<-9- 



6<t^J/V)< — 



= 1.73pb; 

l.Mpb; 
O.Slpb. 



(28a) 

(28b) 
(28c) 



In NRQCD, the conventional J /ip + LH production includes the both the CS and the 
CO contribution. For the CS process e"*"e~ — )■ J/ip + gg, both the NLO QCD corrections jl^ 



22 



and relativistic corrections 
0.4 ~ 0.7 pb at NLO in as and v'i [1 



23| have been calculated. The cross section is found to be 



20 



22[. The NLO QCD corrections to the CO 
contribution have also been obtained [.2l|. If we choose {0\O{^Sl)\0)'^/^ = (3.04 ± 0.35) x 
lO-^GeV^ (0|C»(^P^)|0)-'/'^ = (^9.08 ± L61) x IQ-^GeV^, which are obtained by a global 



fitting of J ftp production data 4^, the cross section of the CO contribution at NLO will 
be about 0.3pb at yU = 2mc, as(/i) = 0.245. Then the total NRQCD prediction for the 
conventional J/ip + LH production will be about 0.7 ~ l.Opb. Unlike the QED background, 
the cut on 6j/^, for example tt /9 < 6j/^ < 87r/9, will only have a minor influence on the 
conventional QCD cross section a^'"^, because both the CS and CO contribution do not 
depend strongly on 6j/^ in a suitable cut condition 

of 6j/^, the cross section of the conventional QCD process can be in the same order as 



the background cross section. Furthermore, the results in 21| shown that CO contribution 
mainly assemble in the kinematic end point region, while the CS contribution is distributed 
in the whole region of < pj/^ < 4.85GeV, so to study the CO contribution, it can be 
further required pj/^ > 3GeV in the measurement. Such a requirement will reduced the 
CS contribution by about 50%, but has little affect on the CO and the QED background 
contribution. In our calculation, we determine the effective vertices of ilj{2S)'~f*, J/ip'j* and 
ip{2S) J / ip + TTTT hj fitting the experimental data and using the i?- value to represent the 
effective vertex of '~f*qq in the calculation of e~^e~ — )■ J /ip{ip{2S)) + qq cross section, this 
indicates that all the possible important higher QCD correction effects to the background 
are included automatically, which makes the uncertainties of our result very small. Based 
on the above analysis, we think further measurement of the J/ip + LH production with a 
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suitable cut condition of 9j/^ and pj/^ will be helpful to understand the role of the CO 
contribution to the J/ip production mechanism in e^e~ annihilation. 

Recently the complete NLOQCD correction to the polarization of J / ip hadroproduction 



were obtained by two groups 



45 



46| . Due to their different ways of fitting the CO matrix 



elements, they got completely different conclusions. After taking into account the feed-down 
contribution of Xcj and ip{2S) states 4^, the authors found that there is no solution to fit 
the pt distribution of the cross section and J/ip polarization measured by CDF collaboration 
simultaneously. Understanding the J/ip production at B-factories can also help to resolve 
the polarization problem of J/ if) hadroproduction. 

In summary, we study the dominant background sources of J/ip + LH production in 
e^e~ annihilation, which include the ISR process e^e'^ — )■ il){2S) + 7 and higher QED 
process e^e~ — )■ ip{2S) + //, where / can be lepton or light quark, as well as the direct 
e+e^ — i- J/ip + qq process with q = u,d,s quark. We find that the cross section of the 
background process is very large in the whole phase space region. If we make a cut on the 
angle between J/ip and e^e~ beam, the cross section of the QED processes will reduced 
largely and become comparable to the cross section of the conventional QCD process. This 
indicates it is possible to measure the cross section of J/ip + LH from conventional QCD 
production at B-factories. 
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